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Consider  solutions  < H(x,e) ,G(x,c) > of  the  von  Karma n equations  for 


the  swirling  flow  between  two  rotating  coaxial  disks 


Ns  also  assume  that  |h(x,c)|  <_  b/c  while  |g(x,c)|  £ B.  This  work  considers 
the  shapes  and  asymptotic  behavior  as  c ♦ 04.  Ns  consider  the  kind  of  limit 


Sanctions  that  are  permissible.  The  only  possible  limits  (intsrior)  for 


Q(s#c)  are  constants.  Zf  that  limit  constant  is  not  zero,  then 
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SIGNIFICANCE  AND  EXPLANATION 


Under  appropriate  conditions  the  steady-state  flow  of  fluid  between  two 
planes  rotating  about  a common  axis  perpendicular  to  then  nay  be 
described  by  two  functions  H(x,c),  G(x,e)  which  satisfy  the  coupled  system 
of  ordinary  differential  equations 

CH1V  ♦ HH”  ' ♦ GG'  - 0 
tC*  ♦KG'  - H'G  - 0 . 

the  quantity  c > 0 is  related  to  the  kinematic  viscosity  and  ^ * is 

usually  called  the  Reynolds  number. 

These  equations  have  received  quite  a bit  of  attention.  First  of  all, 
people  who  are  truly  interested  in  the  phenomena  siodeled  by  these  equations, 
*•9*  fluid  dynanicists,  are  interested  in  this  problem.  However,  as  these 
equations  have  been  studied  by  a variety  of  mathematical  methods,  they  have 


taken  on  an  independent  interest.  The  major  methods  employed  have  been 
(i)  Matched  Asymptotic  Expansions  and  (ii)  Numerical  Computations.  In  both 
approaches  technical  problems  have  appeared.  There  may  be  "turning  points," 
l«e.  points  at  which  H(x,c)  - 0.  Such  points  require  special  and  delicate 
analysis  within  the  theory  of  (i) . As  numerical  problems,  these  equations 
■fa  "stiff"  - precisely  because  c is  small.  The  occurrence  of  "turning 
points”  only  makes  computation  store  difficult. 

for  these  reasons,  these  equations  have  become  "test"  problems  for  methods 
of  "matching  in  the  presence  of  turning  points"  and  "stiff  O.D.E.  solvers." 
■oewver,  when  one  has  "test  problems,"  one  needs  to  know  the  answers. 
Unfortunately  here  the  answers  are  largely  unknown. 

Zn  this  report  we  study  the  asymptotic  behavior  as.  c becomes  small. 


A wealth  of  qualitative  information  is  obtained  which  will  enable  one  to 
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I . Introduction 


In  1921  T.  von  Karman  (SI  developed  tho  similarity  equations  for  incompressible 


axi-syemetric  fluid  flows.  In  1951  G.  K.  Batchelor  [1]  used  the  von  Karman  approach 


to  study  the  fluid  motion  between  two  rotating  planes,  rotating  about  a common  axis 


perpendicular  to  them.  Despite  the  passage  of  time  and  the  work  of  many  people,  this 


problem  is  far  from  being  completely  understood. 


Following  Batchelor,  K.  Stewart  son  [20]  made  a further  study  of  the  problem  and 


disagreed  with  several  of  Batchelor's  basic  conclusions.  In  the  ensuing  years  nany 


people  have  attacked  this  problem.  Numerical  calculations  have  been  carried  out  by 


and  Rogers  [7],  C.  E.  Pearson  (lSl , D.  Greenspan  [3],  D.  Schultx  and  D.  Greenspan 


[19],  L.  0.  Wilson  and  N.  L.  Schryer  [23],  G.  L.  Mellor,  P.  J.  Chappie  and  V.  K.  Stokes 


[13] , r.  D.  Nguyen,  J.  P.  Ribault  and  P.  Florent  [14],  S.  P.  Roberts  and  J.  S.  Shipman 


(XT]  • Formal  matched  asymptotic  expansion  methods  have  been  applied  by  A.  M.  Watts  [22] 


Mb  also  did  numerical  calculations)  K.  K.  Tam  [21],  H.  Rasmussen  [16],  B.  J.  Tsthowsky 


end  W.  L.  Slepun  [12],  On  doubt  ably  many  others  have  also  worked  on  this  problem  and  we 


era  unaware  of  their  efforts. 


Rigorous  mathematical  results  are  a bit  sparse.  There  are  (to  our  knwledqe)  exactly 
three  papers  concerned  with  the  existence  question,  S.  P.  Bastings  (4),  A.  R.  Elcrat  [2] 


and  J.  B.  McLeod  and  S.  V.  Barter  [10  ].  The  first  two  obtained  existence  and  uniqueness 
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result*  for  Hull  valuta  of  tho  Reynolds  number  (R  ■ ) . The  third  concerns  itself 

only  with  the  csee  of  counter -rotating  planes.  An  existence  theorem  i*  obtained  for 
all  c > 0 and  a complete  asymptotic  description  (as  c ■*  Oe)  is  given  for  the  solu- 
tions obtained.  Not  a word  is  said  about  unicity,  h later  paper  by  McLeod  and  Parter  till 
gives  a negative  result  (in  the  limit  as  c ■*  0)  of  the  existence  of  solutions  which 
are  monotone  in  the  angular  velocity.  That  results  contradicts  a conjecture  of  Batchelor . 

The  interplay  between  all  of  these  approach*  has  been  extremely  profitable  and 
interesting.  The  conjectures  and  remarks  of  Batchelor  and  Stewar tson,  the  "shapes* 
obtained  In  numerical  calculations  and  the  general  qualitative  results  of  the  formal 
asymptotic  expansions  have  all  led  to  "target"  questions.  In  their  turn  these  target 
questions  have  been  studied  analytically,  numerically  and  by  formal  expansion  methods. 

Por  example,  the  results  of  (10]  cast  doubt  on  the  calculations  of  ( 3 ] and  a refined 

method  was  proposed  in  (19].  One  of  the  goals  of  (12)  was  to  obtain  - via  formal 

* 

expansion  technique*  - tho  solution!  of  (10). 

hot  ue  now  describe  the  problem.  Let  the  plane!  be  placed  at  x « 0 and  x • 1 
aad  rotate  about  the  x-axia  with  constant  angular  velocities  0Q,  respectively, 
hot  qr»  q^,  q^  denote  the  velocities  in  cylindrical  coordinates  (r,S,x).  Following 
mom  Karmen  ( S ] and  Batchelor  ( 1 } we  make  the  onsets  that  in  function  of  x 

aloms,  l.a.  there  is  a function  H(x)  such  that 

q^  “ -H(x)  . 

Than,  as  a consequence  of  the  steady  state  Navier  stokes  equations  we  find  that 

q,  - | H*  (x) 

and,  there  is  a function  G(x)  such  that 

qg  * | C(X>  . 

Theme  functions  (H(x).G(x))  satisfy  the  ordinary  differantial  equations 

1.1)  *HiV  ♦ HH"'  ♦ GG*  « 0 

1.2)  tC"  4 MG'  - H'G  - 0 


-a- 


1 


shere  c la  tha  kinematic  viaeoalty.  Tha  associated  boundary  conditions  are 
J*J)  H(0)  ■ H (1)  • 0,  (no  penetration) 

1-0  H'(O)  . M*  II)  - 0,  (no  slip) 

l.S)  0(0)  - 2nQ.  o(i)  - an4  . 

However,  our  results  are  independent  of  these  boundary  conditions.  Hence 

they  apply  to  tha  cases  where  one  has  "suction"  or  "blowing"  on  the  planes. 

In  this  work  we  are  concerned  with  the  asymptotic  behavior  of  solutions 

(Ib,^) ,C(x,ca))  as  ♦ 0-»  under  the  basic  hypothesis: 

1.1)  There  is  a constant  B such  that 

|H(x.c  ) I < bvT  . 
n — n 

|C(x,cb)|  <_  B . 

this  hypothesis  has  been  used,  implicitly  or  explicity,  in  many  of  the  studies 
eomeeeted  with  this  problem.  There  are  good  reasons  for  this.  Tor  example,  if  wm  set 

C • ^ . h(C,c)  - H (x,c) , g(C,c)  - G(x,e)  , 
than  equations  (1.1),  (1.2)  became 

(i]4h  * h(A)Si,  • , A , . 0 , 

*•*'  («)  « * h(s]’  ' («)’  ' 0 • 

em  the  larger  interval  (0,  ^ ] . Many  of  the  formal  asymptotic  expansion  studies  [211, 
Oil  ( 22 1 have  used  this  "stretching"  and  then  "matched"  with  the  numerical  results  of 
snd  Lance  [18]  for  the  semi-infinite  region,  i.e.  the  von  Karaan  problem, 
rical  studies  based  on  "shooting"  methods  (111,  [171  have  found  it  convenient 
taaahe  this  change  and  then  actually  compute  < h({,c)  ,g(C,c)  >.  In  fact,  Wilson  and 
Sshryar  [231  who  did  not  use  a shooting  method  also  found  these  variablee  convenient 
for  eslculation.  Moreover,  the  solutions  found  by  McLeod  and  Farter  (121  do,  in  fact, 
satisfy  these  estimates. 

tee  may  integrate  equation  (1.1)  to  obtain 

!••)  tH"*  ♦ HH"  ♦ ~ GJ  - ~ (H’)a  • V(C) 

-3- 


L 


f 


where  |i <c)  is  a constant.  This  constant  is  of  some  independent  interest,  ror 
sssayle.  In  ths  semi-inf initc  problem  i.s.  single  disk  or  von  Karman  problem,  one  sets 

1.10)  v - j Lim[G(x))2  > 0 , 

* x-«— 

and  v is  a known  quantity.  On  the  other  hand,  in  110)  i_t  was  found  that  u(c)  — -t . 

And,  of  course,  in  the  two  disk  problem,  v(c)  is  unknown. 

Xn  section  2 we  set  C ■ 1 and  study  solutions  of  (1.1),  (1.2)  on  large  intervals 

(0,E)  with  E >>  1.  Here  we  discover  several  points  of  interest  for  the  single  disk 

problem.  Xn  section  3 we  return  to  the  finite  interval  and  values  of  c < 1.  Under 

the  assumption  H.l  we  are  able  to  make  the  change  of  variables  (1.7)  and  apply  the 

results  of  section  2.  The  main  results  of  this  section  are  (i)  If  u(c  ) ■*  w as 

n 

* 0+,  then  |i  >_  0,  (ii)  one  My  select  a subsequence  n^  ■*  0 such  that,  there  is 
• constant  ga  and  for  every  4,  0 < 4 < ~ 

1.11)  Nax{|c(x,c  ) - g |i  0 < { < x < 1 - 4)  + 0 as  c ■*  0 . 

x “k  “ ~ "k 


) • 

i results  are  consistent  with  the  suggestions  of  both  Batchelor  (who  emphasised  the 
9m  * 0)  and  Stewartson  (who  emphasised  the  case  gm  m 0) . 

Xa  amotion  4 we  consider  the  case  u > 0.  Xn  this  case  we  find  that  there  is  a 


•Mutant  hw  such  that 


IH(x 

~7l 


H(x.«  ) 

V 


h-  ' 4 


<x<l-«J-0. 


Mich  of  our  work  described  in  this  paper  is  based  on  the  properties  of  the  function 

t)  #(x,c)  - tC'(x,c))a  ♦ (H"(x,c))2  . 

basic  result,  due  to  McLeod  ( 9|,  (10),  is 

aa  <i  The  function  4(x,c)  satisfies  the  differential  equation 
I)  ♦ H4*  - 2c[(G")2  ♦ CH*»  •)*!  , 


. 1 


Md.  tho  function 


1 

•'(X.c)OX p[~  J H(t.C)dt) 

*0 

in  nondecreasing.  Sine*  it  la  also  holoaiorphic.  it  has  at  noat  one  saro.  Thus'  tha 
btluvior  of  ♦(x.c)  is  described  in  one  of  the  following  three  wayst 
(a)  ♦ is  oonotone  decreasing  on  its  interval  of  definition. 

(I)  • is  Monotone  increasing  on  its  interval  of  definition, 

(T)  there  Is  an  interior  point  y such  that  •*  < 0 for  x < y and  •*  > 0 for 

« > V. 


I 


5- 


r 


L 


2.  Bone  Basic  titiwtti 

In  this  aaction  we  *ra  concerned  with  obtaining  aatlnataa  on  functiona  <h(0,g(0> 
ehlch  satisfy  the  differential  equations 

iv 


h ♦ hh1"  ♦ gg*  - 0,  0<  ( < J 

g"  ♦ hg*  - h*g  - 0,  0 < ( < I 


2.1a) 

1.  lb) 

(inn 

a.l)  0 < ) < i < • . 

Mar aowar , theae  functiona  aatlafy  the  a-priori  estimate 
a.3)  |h(()|  « B,  |g(C)|  < B . 

throughout  thia  aaction  tha  lattara  C,  B will  denote  thaaa  conatanta. 

Wa  recall  a basic  aatinata  dua  to  Landau  ( 6 ) . 

II 

2.1»  Lot  f(C)  • C (0,0  and  lat  n > 0 be  a given  poaltive  number.  There  la  a 


a.d) 


stant  C(n.N)  depending  only  on  i\  and  N and  not  on  the  length  K,  such  thati 

1 < J < N - 1 

Itef'l.  i c«n.N>  ||*||_  ♦ • 


# if  n < f I than 


a.S) 


«*•!.  i?  Ilf  II.  • J II 'll.  • 


t See  (61. 

hi  In  aoat  Instances  this  1« 
use  (2.5)  when  n is  large. 


Is  applied  when  h la  snail,  however  we  shall 


2.2«  Let  h((),g(C)  satisfy  (2.1a),  (2.1b)  and  (2.3).  There  are  constants  B^, 


J • 1*2*...  (depending  only  on  B and  not  on  h,g)  such  that 

**«  111^4  IV 


Proof  i Let  i)  ■ . Prom  (2.1a)*  (2.1b)  and  leana  2.1  we  obtain 

||hlv||.  1 blbc(n,4)  ♦ n||hlv||j  ♦ B(BC(n,2)  ♦ nRfll.). 
Ilf  II.  < B(BC(n*4)  a n||hiv||j  ♦ B(BC(n,2)  a nUg"||.)* 


> ■ 

t 1 


Collecting  tens  end  adding  the  Inequalities  we  obtain 


(l  - a»n)l||hlwH-  ♦ ||g"||j  < 2B2lc(n.4)  ♦ c(n.2)i  . 

That  is 

l|hiV|L  ♦ ll«"IL  I 4B3lC(n.4)  ♦ C(n,2)]  . 

Ttuia,  (2.6)  follows  frosi  (2.4)  and  repeated  differentiations  of  the  basic  equations. 

In  the  remainder  of  this  section  we  use  these  estimates  and  leama  4 to  obtain 
•eon  stronger  estimates. 

Let 

a*T>  ♦ («  - Cg'(C)]2  ♦ (h-(5)lJ  • 

!«■■»  4 (with  c - 1)  applies.  Suppose  4'  (()  > 0 on  a "large"  sub  interval 
•f  [0«E] . Since 

0 <4(0  ibJ  ♦ b2  . 

them  4'  (O  must  be  "small"  on  "relatively  large"  sets.  Our  next  result  makes  this 
•tatement  precise.  The  details  of  the  proof  are  left  for  an  appendix. 

V— a 2.3»  Let 

>*••>  - *««(i.IUIL.II*,IL.I|4"!I.}  . Ko  " *0  * W ' 

K1-av,,ii.- 

Let 

*•*)  16K^  « l < E . 

ttta  for  every  Interval  (a,0)  C [o,E]  of  length  L,  i.e., 

2.10)  0 - a - L , 

■set  that 

2.11)  ♦’«)  >0,  l * [o.ei  . 

there  is  a subinterval  (a',0']  c (a,0]  such  that 

2.1»a)  §•-«•»  -i—  lV2 

16*J 
1 1/4 

2.12b)  0 < 4'(C)  < (f) 


7- 


Moroover , on  this  Interval 


2. 12c) 


!♦•«)!  I 0^  ♦ l)  (i)  /e  . 

froot i The  estimates  (2.12a).  (2.12b)  follow  inoicdlatcly  from  Theorem  A of  tho  Appendix 
while  (2.12c)  follows  fro*  (2.S)  applied  to  ♦ ' with 

•-nr. 


Corollary  2.3i  On  this  saae  interval  ta',0']  we  have 
2.13)  (h 


2 3 -l /a 

♦ (9*)  * BL  ♦ (1  ♦ A1)L  < K„L 


1/8  < k2i'1/8 


hoof  i Apply  (2.3),  (2.12b)  and  (2.12c)  to  the  differential  equation  (1.16)  (with  c - 1) . 
. 2.4«  Suppose  (o', S']  C (O.E)  is  a large  interval,  i.e. 

_ a'  > — L_  Ll/2 


2.14) 

on  which 

2. 15) 

l is  so  big  that 


♦ |g"|  < v‘1/16 


2.16)  X, 1/64  < L1/32  < L1/2  < j (B*  - o')  . 


32*0 


there  is  a constant  such  that 


2.17) 

2.  It) 


|h-(C)  | « |g‘  (C)  | < k4l“1/32 

|b'(C)  | l k4L*1/64  . 


,1/32 


Hoof  i Lot  n - L . Applying  (2*5)  to  tho  function  h*  (C)  ve  have 


B*»- II.  i \ t1/3Jllh,,,iL  ♦ 2t'1/32lli»,IL 


if  V ♦ -i 


W£i 

»#  with  n “ X.  we  obtain 

H“,IL  I \ I-l/64||>»-||.  ♦ 2L"1/64||h||.  1 K4L-1/M  . 

A similar  argument  gives  the  result  for  |] g * ||  . 


Further  applications  of  lasna  2.1  give  the  following  additional  estimates. 
JjSSLLi1  Suppoae  (a.fl)  C (0,E)  is  an  interval  on  which 


2.19) 
and 

2.20) 


0 < 4(C)  I KL 


-1/16 


. > I L1'” 


then,  there  is  a constant  M such  that 

2.21a)  It*  CO  | iML~1/32  ' 

2.21b)  |4’(C)|  <ml"1/64  , 

2.21c)  |h,"(C)|  < Ml‘1/64  * 

2.21d)  Jg“(C)  | £HL~1/64  , 

2.21e)  |h’(C>|  <ML-1/64  . 

theorem  2.1»  Suppose  E - +*•.  Then,  either  4 ((}  s 0,  or 
2.22)  4*(()  < 0,  0 < C < • . 

troofi  Suppose  there  is  a point  CQ.  0 < CQ  < • at  which  (2.22)  is  violated.  Then 


4’<0  > 0, 


C0  " C < ♦-  • 


tot  L be  so  large  that  we  nay  apply  Lemma  2.3  in  the  interval  (CQ  ♦ L,  CQ  ♦ 2L] . 
wu  find  an  interval  (a, 6]  C [£Q  ♦ l,Cq  ♦ 2L)  such  that 


|4"(C)|  « KL_1/8 


|b*  "(C)  | < Kl'1/16,  |g"  (C)  | < KL"1/16  . 

If  L Is  sufficiently  large  we  may  also  apply  Leona  2.4  to  discover  that 

*.»)  |4(C)|  £ KL'1/16  . 

*T  *he  neo»re  of  4(0*  this  last  estimate  holds  on  the  entire  interval  IC0«C0  ♦ t] . 
toMever,  since  L is  arbitrary  we  have 

4(0  - 0.  CQ  < C • 


towever,  since  4(C)  is  a holomorphic  function,  4(0  10 


Kenarhi  This  result  is  similar  to  a result  of  McLeod  (111.  However,  tho  proof 
is  quits  different,  as  are  the  hypothesis. 

Theorem  2. 3i  Suppose  C ■ ♦«*.  Then  Lim  g(t)  exists,  call  it  gw.  The  constant  of 

t— 

integration  p is  given  by 

1 3 

» • J • 

►roof'  Choose  a large  number  L and  let 

MO  - -h(3L  - {),  0 < ( < JI, 

g(()  - g(3L  - O.  o < t « 3L 

4(0  - (hH)3  ♦ (q*)J.  0 < t « 3L  . 

Then,  <h((),g(C)>  satisfy  (3.1a),  (3.1b)  and  (3.3).  Moreover,  4(0  satisfies  (with  - 
appropriately  placed)  (1.16).  Finally, 

4*U>  * “4*(3L  - O > 0,  0 < ( < ll  . 

As  in  the  proof  of  theorem  3.1  we  apply  lemma  3.3  on  the  interval  (L,3L).  Applying 

laema  3.4  we  find  that 

0 < 4(()  < ml"1716, 

Applying  lemma  3.5  we  find  that,  as  L •* 

h*  (()  - 0(l"1/M), 
h-(U  - o(l'1/5J). 
h**MO  - o(l'1/64). 

That  is 

hMO  - 0(l’1/64), 
h-(0  - o(l-1/3J), 
h*"(0  - o(l"1/M), 

Inserting  these  estimates  into  (1.9)  gives 
(3.34)  ~ g*(0  * U as  C . 

thus  p > 0.  It  p • 0 then  g_  - 0.  if  p > 0,  |g(())  is  bounded  away  from  sero  for 
( sufficiently  large.  Let  o « sgn  g(C),  C large.  Then  the  theorem  follows  with 

g„  ■ (^5p)o  . 


°<(  <i  . 

0 < ( < L , 

0 < ( < l , 

0 « ( < L . 

L < C < 3L  , 
L < ( < 3L  , 
L < ( < 3L  . 
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3.  The  asymptotic  Drhavior  of  G(x,c) 

Be  turning  to  the  (unctions  <H(x, e) ,G(x, «)>  which  satisfy  (1.1),  (1.3)  on  (0.1) 

vs  consider  their  behavior  ss  * *0*.  Of  course,  we  also  assume  H.l,  l.e.  (1.6a), 

(1.6b).  We  make  the  change  of  variables  (1.7)  and  consider  the  ’‘stretched"  functions 

h(C, t),  g(C, c)  on  the  Interval  (0,  ).  Ne  observe  that 

* n 

(^J'hU.c)  - (*)  2 (£)rH(*,«>  . 

3*l>  r r 

(d«)  ,(t*c)  " <c)2  (jta)  G(x’e)  * 

Usm  3. It  Let  < H(x,c) ,G(x,c) > be  a solution  of  (1.1),  (1.3)  which  satisfies  (1.6a), 
11. 6b).  Let  B^  be  the  constants  of  lemma  3.3.  Let 


C0  “ 3 3 Bl  + M3*B3 


3.3)  )f (c)  | < CQ  . 

Proofi  Using  (3.1)  ws  see  that 


K— s 3.3«  Let  4,  0 < 4 < - be  given.  There  exists  an  c ( A)  >0  and  an  N(4)  > 0 
depending  only  on  6 and  B such  that)  for  0 < c <_  c (A)  and  4<^*Tci.l“4ve  have 

3.3a)  |h«({,c)|  < M{A)c1/64.  |g' (C.c)  | <M(4)t1/M 

*•»»)  |h'(C,c)|  < H((|cl/1,81  |h'  •'«,«)  | ‘mOi171”  . 

hvofi  Let 

4(0  ■ («'((»c)la  ♦ Ih-(C.e) Ja  . 

Them  4(C)  satisfies  (1.16)  with  c • 1 and  leans  6 applies.  Let  Kg,  K^  be  as 
la  lemma  3.3.  Let 

( c(4)  - «a/(33)aKj  . 

then,  if  0 < c < c (4) 

3.4)  16K^  < - L < ^ - E . 
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Ut  y(c)  be  th«  unique  point  at  which  4(0  assumes  its  minimum. 

Cm  li  V(c)  « In  this  case 

♦ •(C)  > o.  i . « < i/T  c « i - 1 < , 

4'(C)  < 0,  J 4 < 4 ( < { . 

Tha  estimate  (3.4)  Implies  that  (2.9)  holds  for  the  two  intervals  JfJ  ' 

Tx  - « 1 - £*1 

I *~Jf~  • — J/f — !•  Hence  we  may  apply  lcima  2.3  and  lemma  2.4  to  obtain  subintervals 
outside  and  * constant  M(4)  so  that  (3.3a)  holds.  Since  4(C)  assumes 

its  ataximum  on  the  boundary  of  any  Interval,  (3.3a)  holds  on  all  of  j . The 


Hence  we  may  apply  lcima  2.3  and  lemma  2.4  to  obtain  subintervals 
-1  and  a constant  M(<)  so  that  (3.3a)  holds.  Since  4(C)  assumes 


estimates  (3.3b)  now  follow  from  lemma  2.2. 


2.  y (*)  r 


For  definiteness  suppose  y(*)  <_ 


Then  we  arque 


Tl-4  1 " j 4 1 

as  above  to  obtain  the  estimate  (3.3a)  on  a subinterval  of  I — ^ ■ I . Then 

■luce  4*  > 0 for  C > Y(»)  we  see  that  (3.3a)  holds  on  ail  of  j • he 

before,  (3.3b)  follows  from  lemma  2.2. 

L— **  3.3»  Let  4,  0 < 4 < j-  be  given.  Let  e £c(4)  and  let  < H(x,t)  ,G(x,c)  > be 
a solution  of  (1.1),  (1.2)  satisfying  (1.6a),  (1.6b).  Then,  for  4 ^ x 1 - 4 we  have 


I jG  (x.c)  - u (t)  | < M(6)[c1/13B  ♦ 


1/128  . „ 1/64  , 1 


♦ ~ H(6)  t ] 


Froof i Me  make  the  change  of  variables  (1.7),  using  (3.1)  and  (3.3a),  (3.3b)  we  obtain  (3.5) 
■ota i While  the  analysis  given  in  this  paper  is  primarily  concerned  with  "limit" 
behavior  and  families  (i.s.  sequences)  of  solutions  ( H(x,c) ,G(x,c) > which  satisfy 
■•1  r (1.6a),  (1.6b)  the  estimate  (3.3)  provides  a "check"  which  may  be  applied  to  any 
calculated  pair  ( H(x,c)  ,G(x,c)  >.  We  simply  must  carry  out  some  messy  computation, 
that  Is,  (i)  find  a B for  (1.6a),  (1.6b)  j (ii)  carefully  follow  the  steps  of  section  2 
and  compute  M(6)i  (iii)  check  (3.5). 

theorem  3.1»  Let  c 4 0>  and  let  (H(x,c  ),C(x,c  )>  be  a corresponding  sequence  of 
* 1 n n n 

solutions  of  (1.1),  (1.2)  which  satlsfios  H.l,  i.e.  (1.6a),  (1.6b).  Suppose 
3*6)  v (c  ) ■*  y as  c ■»  0*  . 

II  II 


* ism--'. 


(Note i In  view  of  lcnuu  3.1  one  con  always  extract  s subsequence  so  that  (3.6)  holds.) 
Ikon 

3.7a)  U 1 0 

andi  for  every  4,  0 * * * 

3.7b)  Max{|c2(x,cn)  - 2u|j  i « x < 1 - <>  ■»  0 as  tn  ♦ 0*  . 

Moreover,  if  u ■ 0 then 

3.9)  Max(|G(x,cn)|>  J<x<l-Jl  + 0 as  CQ  ♦ 0+  . 

If  v » 0 then  there  is  a subsequence  ■*  • and  a square  root,  say  A,  of  2|) 
each  that 

3.9)  Max{  Icfx.c^)  - a|»  (<x<1*«)*0  as  en  ■*•  0+  . 

Zb  fact,  if 

3.10a)  |v(en)  - v|  ♦ M(«)[t^128  + Be*/64  ♦ | M(6)c*/64)  < o < Jq  Z 

3.10b)  ic2(*,tn)  - W **'  «<*<!-«• 

Thus  G(x,e  ) is  of  one  siqn  and 
o 

3.11)  |G(x,c  , - ^ sgn  G(x,e  ) | < . 

n " ~ /2v 

Proof i The  estimates  (3.7b)  follows  imnediately  from  (3.5).  The  inequality  (3.7a) 
fallows  from  (3.7b).  Then  (3.8)  is  apparent.  When  v > 0 and  (3.10a)  holds  we  have 
frcei  (3.5)  and  the  triangle  inequality 

I 2 ^ -®  - W " * 

(3.10b)  follows  at  once.  We  then  have  • 

|c(x,c  ) - &\  • | G (x , c ) ♦ i® 
n n 

aad  (3.11)  follows  at  once.  Thus,  choosing  "signs"  at  a fixed  point,  say  xQ  " j » 
we  obtain  (3.9). 
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4.  The  Asymptotic  Behavior  of  H(x,c)  : u > 0 
Let  4<  0 < 4 « be  given.  Let  t(4)  be  the  value  determined  in  soction  3 
and  suppose  0 < c « c<4).  Let  ( H(x,r ) ,G(x,c)  > be  a solution  of  (1.1).  (1.3)  on  (0,1) 
which  also  satisfies  H.l,  l.e.,  (1.6a),  (1.6b).  We  also  suppose  that  there  is  a 


constant  u > 0 such  that 


V(c)  > u/2  > 0 , 


j C2(x.t)  > 5/4  > 0.  j 4 < X « 1 “a®  * 


The  stain  result  is 


theorem  4.1i  There  are  positive  values  c - t(4,u).  K - K(4),  o - o(u)  where  c depends 
only  on  6,.  p and  B,  K depends  only  on  4 and  B,  while  o(y)  depends  only  on  ii 


and  B,  such  that,  for  4 x <.  1 - 4 


|±Ga(x,c>  - u(c)|  <_  K(4)exp{-oc  l^584) 


|(to)  i * 2 K(4)axp{-oc"1^3M), 

|jjG(x,c)|  < c-1/2K(4)exp{-oc"V384)  . 


-1/384,  , , , 

-oe  ),  r • 1,2,3  , 


Finally,  there  are  constants  a,  b such  that 

4.3a)  | H(x,c)  - a|  < c 1^2K(4)exp{-oc  , 

4.3b)  |G(x,c)  - b|  < t_l/2lC(4)exp{-oe"1/384}  . 

4.3e)  |b|  ♦ /2ti(c)  . 

Of  course  this  theorem  inn ed lately  implies  certain  limit  theorems  for  subsequences 
of  solutions. 

The  proof  is  relatively  straight-forward  and  follows  the  general  pattern  McLeod's 
work  in  ( 8 ] . Unfortunately  there  are  many  details  to  check  out.  We  outline  our 
approach. 

I5SJLA'  **•  the  change  of  variables  (1.7)  and  consider  the  functions  h(C,c),  g(C.c) 

on  an  "interior"  interval  (a, 8]  which  satisfies 

4.4a)  4//T ia  < I < (1  - t)//t  * 

4.4b)  I - a - L(c)  - c”1/384  . 


ii 

. . e 


SSSEJL*  **•  ,ind  *n  equivalent  integral  equation. 

|tsE_3.  Ne  prove  local  unlqueneaa  for  solutions  of  the  inteqral  equation. 

•teg_4i  Ne  prove  that  the  desired  solution  of  the  integral  equation  can  be  obtained 
via  Heard  Iteration. 


•£eg_5»  *•  that  the  limit  of  the  Picard  iterates  satisfies  the  appropriate  estimates. 

JtegJSt  Ne  return  to  the  original  variables  x.  H(x,c)  .C(x,c) . 


Ne  can  imagine  step  1 has  been  done. 
Use.  2*  *h«  Integral  equation.  Let 


4.Sa) 

4.9b) 


hQ  - h(a.c) , 

«X(C)  - g(C.c)  - «Q, 
«,<C>  - g'(C.c). 
o5(C)  - h-(C.c). 


g0  - C(o,c)  , 

u2({)  - h(*,c)  - hQ 
u4(C)  - h*(C,e)  . 
u6(«)  - h' ' * (C,c)  . 


Let 

4.4a) 


4.4b) 


0 - «WVVW 


0 0 1 


0 0 0 


0 0 0 1 0 

0 0 -h  g 0 

0 *0 

0 0 0 0 1 


0 

0 

0 


0 0 0 0 0 1 


0 


0 


■*0 


0 


0 


4.4e)  b ■ b(0)  ■ (0,0, u.u  * u.u  ,0,0,  - U u -uu)*. 

• ^ 23  1 3 26 

fba  equations  (2.1a),  (2.1b)  now  take  the  form 


«•*>  dC  " + b(0)  * 

b direct  calculation  shows  that  the  eigenvalues  of  A are  the  roots  of 

*•*>  »*t«J  ♦ ♦ h0)J)  - 0 . 
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kssentially  this  hm  eigenvalue  problca  arises  in  ( 6 ) and  w can  easily  check  the 
following  formula*  for  the  eigenvalues  1^,  k - 1, 2,3,4. 5, 6.  Let 

•i  " *a  ho  4 VS  ((hJ  4 4 ho)1/J  4 0 • 

"a  " "a  ho  - VS  (<hJ  4 4 hJ)1/J  « 0 • 

» - a 7t  {,hJ  4 - hJ)VJ  * 


Thee,  the  eigenvalues  of  A are 
4.9a) 


* Xj  ■ 0 , 


4.9b) 

4.9c) 


X,  • Pj  ♦ it. 


x5  - fi3  ♦ it. 


PA  - it  . 


Xj  • Pj  “ it  . 


It  la  easy  to  see  that  P4  > 0 and  * 0 provided  that  gQ  * 0.  However, 
(4.1b)  gives 


»*  > «o  I v/i 


4.10e) 

ted,  for  all  « , we  have 
A 

4.iob)  la^l  < » . 

Mm  a simple  compactness  argument  show*  that  there  is  a constant  p > 0 such  that 
4.10e)  *]1*P<0<P1P^  • 

I*t  us  diagonal lse  the  aatrix  A.  Ns  construct  the  aatrix  of  eigenvectors.  Let 


4.11a) 

aad  let 


■<X)  - 


!♦ 


4.11b) 


T • 


0 

1 

0 

0 


>(!,) 


»(X4) 


■ (X5) 


i(Xt) 


XjbdjJ  l4s»(X4)  X$a(l5)  X%a(X4) 


0 


0 


Osing  the  fact  that 


4.11e) 


*<V 


Xk*h0 


Xk<Xk  + V 


3 < k < 6 


om  can  easily  verify  that  the  colunns  of  T are  eigenvectors  of  A and 

-1 


T AT  » diagonal (0,0, Ij.l^.Xj.l^)  5 A . 


4.12) 


**- 


4. lid) 

Moreover 

“l  0 x x * 

0 1 x x x 

0 0 x x x 

0 0 x x x 

0 0 x x x 

0 0 x x x 

*x"  Barks  an  ele—nt  we  do  not  need  to  compute.  (See  ( 8 ] where  an  analogous 
calculation  is  carried  out.) 
bet 

4.13a)  O - TV,  V - T_10  . 

Then  (4.7)  takes  the  form 

4.13b)  fj[  - AV  ♦ d(V)  , 


4.13c)  «(V)  - T_1b(TV)  . 

Thus,  we  have  essentially  found  our  integral  equation. 
lassM  4.1:  Let  V(C)  be  any  solution  of  (4.13b).  Then 


4.14a)  Wj(«  - VjCa)  ♦ / djWtndt.  J « 1,2  , 


!,«-*> 

!,(€-«  X -X,(t-8) 

4.14b) 

Vj(t)  - ^(Wa  3 

♦ e3  J e 3 d.(V(t))dt, 

ft  3 

j - 3.4  , 

v,cc-«> 

*,(C-a)  t -X.(t-ol 

4.14c) 

Vj(C)  - Vj(a)e  3 

♦ a 3 /a  3 d.(V(t))dt, 

• 3 

j - 5,6  . 

kroof » Integrate  (4.13b). 
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■■■■■■ 


K«P  3i  bocal  uniqueness. 

**  sasential  that  wc  distinguish  between  the  components  v^(()  of  V({).  Let 


4.1Sa) 

HvtOll 

- max(|v  (C)  |»  j - 

1 3 

1,2, 3, 4, 5, 6. ) 

4.13b) 

h(V(t)) 

- »a*{ |v . (()  1 1 J - 

j 3 

3, 4, 5. 4.) 

4.13c) 

MCVtC)) 

- max(|vl(t)  |,  |v2(C)  | 

D • 

Usma  4.2: 

There  are  positive 

constants  s^.a^.a^ 

such  that)  if 

4-vectors, 

4.14) 

| ||d(V)  - d(Y)  ||  £ a^M(V  - Y)  • N(V)  ♦ 

l aa*(Y)N(V  - Y)  ♦ SjNCV  - Y)  (N(Y)  ♦ N(V)J  . 

constants  *re  uniformly  bounded. 

[i  The  coefficients  of  T are  bounded  functions  of  g^hg.  A compactness  argument 
Amm  they  are  uniformly  bounded.  The  form  of  T shows  that!  for  k • 3,4,5, 6 and 
1 - 1.2 

" WVWV  * ®jkCVVW 

tJfc  is  Unear  and  homogeneous  while  is  quadratic  and  homogeneous.  Thus, 

t4.lt)  follows  from  the  form  of  b(U). 

faff*  4-3.*  L*t  **(4)  be  the  constant  of  lemma  3.2.  Let  be  a uniform  bound 

h'l.  u. 

4.17a) 


on 


4.17h) 


4.17e) 


that 


that 


4.1«) 


rlM(4)t1/19a  - 0 . 
100  • - 0 , 


tax  ♦ aa  ♦ 2as)L(c)e  £ j . 


*•4  (8) « Vj(a),  *j(o)  be  specified  so  that 

i{|w1(a)|.|»2tn)|,lvJ(fl)|,|v4(6)|,|vsto)|,|v6(o)|}  <4  . 
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*•••  there  ia  at  aost  one  solution  V(C)  of  the  integral  equation  (4.14a), 
(4.14c)  with  these  boundary  values  and  satisfying 


4. It) 


llvcci  ||  < a . 


Proof  i 
4.30a) 
4.20b) 
fros  lei 
4.21a) 


Suppose  Y((),W(()  are  two  such  solutions.  Let 

D £ Mx(||d(Y(())  - d(W(C))||,  «<(<(} 

B - aax(|(Y(C)  -W(C)||»  «<?<(}. 

■a  4.2  and  (4.19)  we  have 

D - (*1  * a2  ♦ 2a3)E9  . 


(4.14b), 


Ttom  the  integral  equation  we  see  that 
4*21b>  B < L(c)0  < L(c)  (aA 


+ a2  + 2a3)E6  . 


thus,  either  B - 0 or 


1 < L(C)  (a3  + a2  ♦ 2»3> 0 

idlieh  contradicts  (4.18). 
top  4:  the  Iteration. 

***  *1^'  vj(aJ*  v4(B),  v5(o).  vg(a)  be  determined  from  h(£,c),  g(£,c) 

r^*  transformstions  (4.5b),  (4.13a).  We  seek  to  recover  the  appropriate  V{£)  via 
Picard  iterations.  That  is,  let  V°({)  - 0.  Assuming  that  V1^)  has  been  computed 
m determine  Vr+1(()  from  the  equations 


4.22a) 

o$*C« 

- v (o)  + / d (Vr(t))dt,  j - 1,2  , 

a 3 

4.22b) 

▼J+1U> 

x.te-B)  s -i.(t -e) 

-▼.(8)0  +e3  f e 3 d.flr'uMdt, 

8 3 

3 - 3,4 

4.22c) 

^l(« 

MC-o)  l.(C-o)  C -X.(t-o) 

” vj(°)e  / • 3 d.tv'itMdt, 

a 3 

3 - 5,6 

Xeama  4 

^4i  Let 

v'co  be  computed  as  above.  Assume  that 

4(ax  ♦ a2  ♦ 4a3)8V4  < l . 

4.23a) 

4.23b) 

2L(c)8V4  < i , 

Uti 


4.244) 


4.24b) 


4.24e) 


4.244) 


MtV^C)  - Vr‘l(C))  < , 

2 

r r-1  4 T 

*(Vr(C)  - vr  l(C),  < |e  ♦ • ] . 


«(Vr(C))  26  . 


mvrU)) 


i Ml*  ♦ e 2 J . 


H**3"**!  the  functions  V<r) ({)  converge  uniformly  to  4 function  V<{)  which  satisfies 


integral  equation  (4.14a),  (4.14b),  (4.14c)  and 


4*as*>  M<V(t>)  < 26  , 

a Mki  rb|<C-6)  P.U-s)! 

4-Wb>  H(V(C))  < 26  [e  1 ♦ e 2 J . 

2*25*1  Ha  observe  that  (3.2)  together  with  the  choice  of  L(c)  implies  that  the 
notation  V(C)  determined  by  < h((,c)  ,g (C,c)  > satisfies 

4'M)  ||v(C)  ||  <6  . 

tborsfore,  a-fortiori,  the  boundary  conditions  satisfy  the  same  estimate.  Thus,  (4.24a>. 
l4**4^  *r*  eetisfied  for  r ■ 1.  We  proceed  by  induction.  Assume  that  (4.24a), 

M.24b)  are  satisfied  for  r - 1,2 rQ.  Wien.  (4.24c).  (4.24d)  are  also  satisfied 

<or  r " 1*2 rQ.  Applying  lemma  4.2  we  have 

f iion-So,  - d(vv‘<(nii  < -ai.  Ui . ,2 . <v[.v'-#l  4 


4.25b) 


#7/4  rp  (H)  P,(C-oH 


titution  into  (4.22a)  gives 


J • 1,2 


' using  (4.23b)  we  have  (4.24a)  with  r - rQ  * 1.  Substitution  of  (4.27)  into 


(4.22b)  gives 


t 


< 

I 

> 


-L 


"*1 


next  theorem  ia  an  immediate  consequence  of  theorem  4.1. 

Theorem  4.2i  Let  ( H(x.tn)  .Clx.r.^  > be  a sequonce  of  solutions  of  (1.1),  (1.2)  which 

satisfy  H.l.  Suppose  these  are  constants  C.  qm  with  » 2u  such  that 

*•31)  V(c  ) ■*  i!  , 

n 


I,  for  every  4,  0 < 4 < — we  have 


4.32) 


{ |G(x,c  ) - ? |i  « < x < 1 - 4)  ♦ 0 as  c •*()♦. 
n • — — n 


Than  there  is  a subsequence  c -*0+  and  a constant  h so  that 

\ 0 

4.33) 


{ |~^F — H(x,e  ) - h |;  4 <_  x ^ 1 - 4}  ■*  0 as 


c -*•(>♦ 
*k 


Moreover,  if 

C ■+  *rt 

"k 

4.34a) 

than 

4.34b) 

{ I jjjf — H(x,e  > - h(x)  | ; 4<_x<_l-4}-*0 


as  c ■*  0 , 


h(x)  t const  . 


•i 
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Appendix 


In  this  appendix  we  prove  the  following  very  plausible  result:  if  4(C)  is  a 
oth  function  defined  on  a very  large  interval  and  if  4(C)  is  both  positive  and 
nonotone , then  there  are  relatively  large  intervals  on  which  ♦'(C)  is  small.  Unfor- 
tunately  the  complete  proof  is  technically  complicated. 

Theorem  A:  Let  4(C)  satisfy 


A.l) 
A. 2) 
Let 
A.3) 


o < 4(0.  o < e < l , 

♦ '(t)  1 o,  0 < C < L . 

*0  -M*{i.n*ii-.iU'iL#ii4-u-} 


and  suppose  that 

A.4) 

Then,  there  is  a subinterval  [o'.B'l  C [o,L]  such  that 


A.S) 


-a'  * 


16KJ  + 1 


,1/2 


1/4 


and 

A.6J  0 < ♦•  (t)  < [±J 

■a  require  a basic  estimate  based  on  the  mean  value  theorem. 
*•!'  t * C*(0,Ll.  Suppose 


A.7) 

and 

A.$) 

let 
A. 9) 


< M , 


A. 10) 


|H  vUa>o  • 


b - *in{-,  CQ,  1-  CQ>  . 


n (t>  «o'b^i'0+b 
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Wo  construct  a soguonce  of  points 


8Ko 

tissL  of  Theorem  Ai  Lot  I— I " ® i ~£~  • 

with  tho  following  properties 

♦*<v  - ** 


♦ M*l>  - i «lf 

♦ Mx2)  - 4, 

♦ Mx2Jn,  -I4, 

♦ Mx2J+2)  - 4. 


♦ ’(X) 

> 

4 

2 ’ 

for 

♦ ’(X) 

< 

4. 

for 

♦ •(x) 

> 

4 

2* 

for 

♦ Mx> 

< 

4. 

for 

x0  i X < Xj  . 
*i  i x < x2  , 


X2j  ~ X < x2j*l  ' 
X2j+1  - X - X2j+2 


*a  accomplish  this  we  proceed  as  follows.  If  4 * (0)  >_  4 then  xQ  - 0,  if  not  xQ 

la  the  first  point  at  which  ♦(x^)  - 4.  Let  x^  be  the  first  point  larger  than  x^ 

•uch  that  4(x^)  ■ j 4 and  so  on.  If  x^  _>  •j  L then  the  theorem  is  true.  Aasuste 
< ^ L.  By  Leema  1 the  number  of  intervals  Is  finite.  Let  N be  the  last  index. 

Thee  Xjj  ^ L.  If  K even,  then  ♦*(x)  21 J for  x^  < x < L.  Thus 

*D  i *(V  -2  (L  " V * 

That  la 

“•-*»>  if  "ji;  • 

It  N la  odd,  then  ♦*{x)  < 4 for  x^  £ x <_  L.  Thus,  we  can  assume  |l  - x^J  < ~ . 

Therefore 

Let  A be  the  number  of  Interval  <x2j*x2j+1>  “ ©n  which  ♦,(x)  > ~ Me  first 
seek  a bound  on  A.  By  Lemma  1.1 

t*aiu  " X2J I - 2K^  * 
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Thu* 


it  * aj*i 

o i *3) 


«d 

R.ll)  R < 4K*/«2  • 

Similarly,  the  total  length  L*  of  these  intervals  satisfies 

a *aj*i 

*•'  f<W  ♦,“tIK0 

j x2j 


R.ia)  t*  i 2k0/«  < \ . 

The  number  of  intervals  “ on  whlch  ♦ * Cx)  < d - is  (R  ± 1)  and  their 

total  length  L“  satisfies 

R.13)  *•"  > | - *•'  > | • 

Thus 

. , L td2 

j aJ  aj‘i  4<R  + 1)  “ 4(4kJ  ♦ «a, 

which  proves  the  theorem. 


» 
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